Libraries

#every census tract AND every retailer sales
acs_retailers <- readRDS("../../../data/processed/03_join_retailers/acs_retailers_w_demo.rds") %>% 
  select(state:total_sales, joint_latitude, joint_longitude)

#tract-level sales info
acs_tract_sales <- readRDS("../../../data/processed/03_join_retailers/acs_tract_sales.rds")

states_with_sales = c("Arizona", "Louisiana", "Michigan", "New York", "South Carolina", "Texas")

pull geometries

Binding geoms

Making map for all retailers but with clustered concentration

retailer_demo_map